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Abstract 

The application of the max-algebra to describe queueing systems 
by both hnear scalar and vector equations is discussed. It is shown 
that these equations may be handled using ordinary algebraic manip- 
ulations. Examples of solving the equations representing the G/G/1 
queue and queues in tandem are also presented. 

1 Introduction 

Max-algebra [U [2] is the system (M U {e}, ©, (g) , where 

e = — oo, x © 2/ = max(3;, y), x®y = x + y Vx, y E M. 

It has the following properties which can be easily verified 

Vx,y, zGM X © (y © z) = (x © y) © z, x © y = y © x, 
X (8) (y (8) z) = (x (8) y) (8) z, x ^ y = y ^ x, 
X (8 (y © z) = (x (8 y) © (x © z), 

C X , X KX) X X ^ 

X © e = X, where e = 0. 

In the max-algebra these properties allow ordinary algebraic manipulation 
of linear expressions to be performed under the usual conventions regarding 
brackets and precedence of © over © . Moreover, the scalar max-algebra is 
extended to the max-algebra of vectors in the regular way. To emphasize 
parallels between conventional linear algebra and the max-algebra, similar 
notations are used for the iterated operations © and © 

n n 

^ ^ Xi = Xi ® • • • (B Xn, — (8) • • • Xn- 

i=l 1=1 
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The max-algebra theory is currently under investigation. There are a 
number of classical algebraic results reformulated and proved in this alge- 
bra. Specifically, the eigenvalue problem has been solved, and analogues of 
Cramer's rule and the Cayley-Hamilton theorem have been found (see survey 
papers [UE])- Moreover, as a research tool in studying practical problems, 
the max-algebra finds expanding applications in many fields of operations re- 
search and optimization, including the analysis and performance evaluation 
of discrete event dynamic systems [H [3l |2] . 

Although max-algebra models were successfully applied to investigate 
certain classes of discrete event dynamic systems, the models of queues have 
received little or no attention. The purpose of this paper is to show how 
queueing systems may be described using the max-algebra approach by lin- 
ear algebraic equations. As illustrations of handling these algebraic models, 
the solutions of the equations representing the G/G/1 queue and queues in 
tandem are presented. 

2 A Linear Algebraic Model for the G/G/1 Queue 

We start with the linear max-algebra representation of the G/G/1 queue 
which provides the basis for more complicated models of queueing systems. 
In the analysis of queueing systems, it is common to apply recursive equa- 
tions to describe their dynamics analytically. Such equations are normally 
written in terms of recursions for the arrival and departure times of cus- 
tomers, and involve the operations of maximum and addition [3llll|5]. 

To set up the equations that represent the G/G/1 queue in the ordinary 
way, consider a single server queue with infinite buffer capacity. Once a 
customer arrives into the system, he occupies the server provided that it is 
free. If the customer finds the server busy, he is placed into the buffer and 
has to wait until the service of all his predecessors completes. 

Denote the inter arrival time between the A;th customer and his prede- 
cessor by Ofc , and the service time of the A:th customer by . Furthermore, 
let A{k) and D{k) be the arrival and departure times of the kih. customer, 
respectively. As is customary, we assume that Ofc > and > are 
given parameters, whereas A{k) and D{k) are unknown variables. With 
the conditions that the queue starts operating at time zero and it is free at 
the initial time, one can readily represent the system dynamics by the set 
of equations O [U [5] 



Let us now replace the usual operation symbols by those of the max- 
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algebra and rewrite the equations in their equivalent form as 

A{k) = ak®A{k-l), (1) 
D{k) = Tk®{A{k)®D{k-l)). (2) 

Under the properties of the operation © and ® these equations could be 
handled much as if they were ordinary linear equations in the conventional 
algebra. Specifically, applying a usual technique to solve the equations for 
the unknown variables A{k) and D{k), we get 

k 

A{k) = ai (8) • • • (g) Ofc, D{k) = ai • • • (g (g) (g) • • • (g) r^. (3) 

1=1 

2.1 The Matrix Representation 

To produce a matrix representation for the G/G/1 queue let us first de- 
fine the vector D{k) = {DQ{k), Di{k))'^ with components Do{k) = A{k) , 
Di{k) = D{k), and replace the symbols ak and by rofc and Tik, respec- 
tively, /c = 1,2,.... It is convenient to preassign -Do(O) = -Di(O) = e, and 
D{){k) = Di{k) = e for all A; < 0. With the new notations, the equations 
([T][2D may be rewritten as 

Do{k) = Tofc0Do(A:-l), (4) 
Di{k) = nk^{Do{k)(BDi{k-l)). (5) 

Substitution of dH) into (0) and the implementation of distributivity of (g 
over © give 

Doik) = Tofc (gZ)o(A: - 1), 

Di{k) = Tik(^Tok®Do{k-l)(BTik<^Di{k-l). 
We may now represent the model in matrix notations by the equation 

D{k) = Tk0 D{k-1), (6) 
where the transition matrix is defined as 

V nk ® Tofc Tifc J 

3 Linear Models of G/G/1 Queues in Tandem 

In this section we extend the max-algebra linear models to cover systems of 
G/G/1 queues operating in tandem. As the basic system of this type, we 
first consider a series of n queues with infinite buffers. Each customer that 
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arrives into this system is initially placed in the buffer at the 1st server and 
then has to pass through all the queues consecutively. Upon the completion 
of his service at server i, the customer is instantaneously transferred to 
queue i + I, i = l,...,n — 1. The customer leaves the system after his 
service completion at the nth server. 

For the tandem queueing system the equations ([1115]) can be easily gen- 
eralized as 

Do{k) = Tok0Do{k-l), (7) 
Di{k) = nk0{Di-i{k)(BDi{k-l)), i = l,...,n. (8) 

where Di{k) and Tik denote the departure time and the service time of feth 
customer at server i, respectively. 

Let D{k) = {DQ{k) , . . . , Dn{k))'^ be the vector of the fcth departure 
times in the system. Similarly as in the case of the G/G/1 queue, we may 
write the vector equation representing the tandem queueing system in the 
form Q with the lower triangular transition matrix 

/ Tok e e ... e \ 

Tik (8) Tofe Tifc e ... e 



Tn-lk ® ■ ■ ■ ® TQk r„_ifc (g) • • • (g) Tifc Tn-lk ® ■ ■ ■ ® T2k £ 
V Tnk®---®TQk Tnk®---®Tik Tnk ® ' ' ' ® T2k ■■■ r„fc / 

Furthermore, we may find the solution of the set of recursive equations 
(IJllSD as an extension of With usual algebraic manipulations, it can be 
arrived at ([Ij) 

(h i-2 \ 

n® ^ n» ® • • • ® n® > ^ = 1-2, . . . 
j=l j=ii j=i„ j 

3.1 Tandem Queues with Finite Buffers 

Suppose now that the buffers of servers in the tandem system described 
above have finite capacity. The feature of queueing systems with limited 
buffers is that their servers may be blocked according to one of the block- 
ing rules |3]. In this paper we restrict our consideration to manufacturing 
blocking which is most commonly encountered in practice. Under this type 
of blocking, if the ith server upon completion of a service sees the buffer of 
the (i + l)st server full, it cannot be unoccupied and has to be busy until 
the {i + l)st server completes its current service to provide a free space in 
its buffer. 

Consider a queueing system with n servers in tandem, and assume the 
buffer at the ith server, i = 2, . . . , n, to be of the capacity bi, < ftj < oo. 
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We suppose that the buffer of the 1st server, as the input buffer of the 
system, is infinite. Since the customers leave the system upon their service 
completion at the nth server, this server cannot be blocked. 

It is not difficult to understand that the kth completion time at the ith 
server, i = l,...,n — 1, can be represented in usual form by the recursive 
equation [3l H] 

Di{k) = max(max(A-i(^), - 1)) +Tifc, - - !))• 

Using max~algebra notations, the complete set of linear equations describing 
the finite buffers tandem queueing system with manufacturing blocking is 
written as 

Do{k) = Tok(^Do{k-l), (9) 
Di{k) = Tik®{Di^i{k)®Di{k-l))®Di+i{k-hi+i-l), 

i = l,...,n-l, (10) 

Dn{k) = Tnk(S){Dn^l{k)(S)Dn{k-l)). (11) 

Although it is evident that we are dealing here with a linear model once 
again, handling the model in its general form ([Slllip requires rather cumber- 
some algebraic manipulations. Therefore, let us consider more thoroughly 
a simple example of a system with n = 2 , 62 = . The equations pillip are 
reduced to 

Doik) = Tok^Do{k-l), (12) 
Di{k) = Tik0{Do{k)(BDi{k-l))®D2{k-l), (13) 
D2ik) = T2k^{Di{k)(BD2{k-l)). (14) 

Going to matrix notations, we arrive at the linear equation 

D{k) =fk0D{k-l), 

with 

fk = 




The above representation of the transition matrix for the system with two 
servers is easily extended to the case of the system with n servers and 6j = , 
z = 2, . . . , n 



Tk 



( TQk £ e ... e \ 

Tlk ® Tok Tlk e e 

Tn~lk (8) ■ ■ ■ (8) Tofc Tn-lk (8 • • • (8 Tifc T„__ife • • • (g) r2fc 6 

r„A. ® • • • (g) Tok T„fc (g) • • • (g) Tlk Tnk®---®T2k ■■■ Tnk / 
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Note that the matrices and differ only in elements of the upper 
diagonal adjacent to the main diagonal. In these elements become equal 
to e, excluding the one of row which remains equaled e. 

Now we return to the example so as to present the solution of the recur- 
sive equations p2m4p . After usual algebraic manipulations one can obtain 

k I i k 
D2{k) = T2k®Di{k). 

3.2 Closed Systems of G/G/l Queues 

Consider a closed tandem system of n queues with infinite buffers. We 
assume that the customers after their service completion at the nth server 
return to the 1st server for a new cycle of service. There are the finite 
number of customers circulating through the system, at the initial time all 
the customers are placed in the buffer of the 1st server. 

Let us denote the number of customers in the system by c. With the 
condition Dn{k) = e for all < 0, we may modify dZHS]) to write the set of 
equations for the closed system without specifying DQ^k) , in the form 

Diik) = nk(g){Dn{k-c)eDi{k-l)), 

Di{k) = Tik®{Di^i{k)®Di{k-l)), i = 2,...,n. 

To produce a matrix representation, we now define the vector of depar- 
ture times as D (k) = {Di{k), . . . , Dn{k))^ . The vector equation associated 
with the closed tandem queueing system with c customers is represented as 



where 



D{k) = Rk(S) D{k - 1) 



5fc(g) D{k-c), 



Rk 



Sk 



\ Tnk 

■ 

e . 
\^ ■ 



T2k 



Tlk 

' Tlk 

' T Ik Tnk ' 



£ 
T2k 

T2k 



Tnk ) 



e Tnk 



Tlk \ 

' Tlk 

'Tlk j 



In conclusion, consider an example of a closed system with n - 
We have the following linear equations representing this system 

Di{k) = Tik®{Di{k-l)®D2{k-2)), 

D2{k) = T2k(8){Di{k)eD2{k-l)). 



2, c = 2. 
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Traditional methods of solving linear recursions give the solution 

k-2 



Dl{k) = Til (g) JJ^(Tii+i e T2i) ® Tifc, 
i=l 
k-1 

D2{k) = Til (g JJ^(Tii+i e T2i) ^ T2k- 



1=1 
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